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Abstract 

We consider a simple model for an SNS Josephson junction in which the "normal metal" is 
a section of a filling-factor u = 2 integer quantum-Hall edge. We provide analytic expressions 
for the current /phase relations to all orders in the coupling between the superconductor and the 
quantum Hall edge modes, and for all temperatures. Our conclusions are consistent with the earlier 
perturbative study by Ma and Zyuzin [Europhysics Letters 21 941-945 (1993)]: The Josephson 
current is independent of the distance between the superconducting leads, and the upper bound 
on the maximum Josephson current is inversely proportional to the perimeter of the Hall device. 

PACS numbers: 74.45. +c, 74.50.+r , 73.43.Jn 
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I. INTRODUCTION 



The zero- voltage Josephson current in a supercondictor/normal-metal/superconductor 
(SNS) junction jl| arises from Andreev scattering ^ at the SN and NS interfaces. In the 
ideal case, an electron incident on one superconductor from the normal metal will be reflected 
back into the normal metal as a hole, and this hole, on striking the second superconductor, 
will be reflected back towards the first superconductor as an electron. When the relative 
phase of the order parameters is such that constructive interference occurs, the back-and- 
forth process continues ad infinitum and transfers two electrons from superconductor to 
superconductor in each cycle [SKZl]. A round trip takes time 2W/vf, where vp is the Fermi 
velocity and W the separation between the superconductors. The current will therefore be 
evp/W for each open transverse channel. In practice, the probability of Andreev reflection 
is less than unity |8|,l9| and the motion in the metal may be diffusive, but evp/W per channel 
remains an upper bound on the critical current. 

An interesting question arises as to what happens when the "normal" metal consists of 
the chiral fermions at the edge of a quantum Hall (QH) bar 10|. In this case the holes 
move in the same direction as the electrons, so conventional Andreev retro-reflection is 
impossible. A two-electron charge transfer requires a (phase coherent) passage around the 
entire perimeter of the Hall bar, and this lengthy excursion suggests that the small 'W of 



the conventional junction be replaced by t 
perturbative study of a S-QH-S system in 



le much larger perimeter L of the Hall bar. A 



Uj supports this conclusion and estimates that 



the maximum Josephson current will be very small — in the order of 1 nA for mm scale 
devices. In view of ongoing experiments on quantum-Hall Josephson junctions, however, it 
seems worth revisiting the problem to see if devices might be engineered to provide larger 
critical currents. 

In this paper we introduce a model of an S-QH-S junction that is simple enough that 
it can be studied non-perturbatively. We obtain analytic expressions for the Josephson 
current /phase relation to all orders in the S-QH coupling, and at all temperatures. Despite 
our greater control over the model, the key conclusions of the perturbative studies in ll|] (see 
also [12]) are unchanged: at filling fraction z/ = 2 an upper bound for the critical Josephson 
current is given by 2evd/L where Vd is the edge-mode drift velocity and L is the perimeter 
of the Hall device. Further, the temperature scale at which the Jospehson current is washed 
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out by thermal effects is set by the edge-mode level spacing En+i — En = 27ihvd/L. Thus, 
if we wish to see Josephson-j unction physics in quantum Hall devices, we should construct 
the junctions by coupling superconducting probes to meso-scale Hall-dots. 

In section two we introduce the model and solve the associated Bogoliubov-de Gennes 
equation. In section three we introduce an analytic regularization scheme to handle the 
otherwise ill defined sums that appear in the current /phase relation. In section four we 
demonstrate that our regularization scheme is consistent with conventional perturbation 
theory at both zero and non-zero temperature. We finish with a brief discussion of effects 
that we have not taken into account, and that may or may not be significant. 

II. THE MODEL 

We consider a z/ = 2 quantum-Hall edge (two spins therefore) in interaction with super- 
conducting (SC) leads (figured]). We model the system by a linear-dispersion edge- mode 
hamiltonian 



Here is the edge-mode drift velocity that is proportional to the gradient of the confining 
potential. The terms with A(a;) are non-zero only where the edge state lies under the super- 
conducting leads. They account for the Andreev coupling arising from the two-dimensional 
electron gas (2DEG) wavefunctions reaching up to touch the superconductor as they drift 
under the electrodes. (See figure El) In contrast to the usual proximity effect, the topological 
protection of the QH edge modes means that this interaction cannot open a gap — but it 
may, for example, convert a charge-(e) right-going spin-up electron into a charge-(— e) right- 
going spin-up hole, and in the process transfer a spin-singlet pair of charge- (e) electrons 
from the Hall bar to the superconductor where they merge with the S-wave condensate. 

We have not included Zeeman-energy term to spilt the energy between the spin up and 
spin down edge modes. Such a term adds only a multiple of the identity matrix to the BdG 
operator, and so has no effect on the subsequent analysis. Further, we assume that the 
energy scales of relevance are smaller than the energy gap of of the superconducting leads. 
We therefore regard the parameters |A| as being externally imposed, and not to depend the 
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FIG. 1: A Hall bar with superconducting probes passing a current I through the edge modes. The 
circled numbers label the regions (1) "outside the leads," and (2) "between the leads." 
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FIG. 2: The wavefunction for an electron in a 2DEG is confined in the vertical direction, but 
there is some amplitude for the vertically oscillating electron to touch the superconductor. As a 
slowly-drifting Landau-level wavepacket passes under the superconducting lead, there will be many 
opportunities for Andreev reEection to turn the electron into a hole. 

energy of the Hall-bar electrons, or on the temperature. 
We can rewrite H in the BdG form 

' -ivd{d^ - ieA) |A(a;)|e^^(^) 



H 



dx < ( ipj^ ) 



|A(x)|e-*''(^) -iVd{d^ + ieA) 



+ const. 



(2) 



Here we have used an integration by parts together with the anticommutation property of 
the Fermi fields to write 



{tpli-ivdid:^ - ieA))ilj^}dx = / {ipi{-Wd{dx + ieA))iljl}dx + const. 



(3) 



This rewriting is essentially a charge-conjugation transformation that makes manifest the 
particle-hole symmetry of the linearized edge spectrum. In particular, it reveals that the 
charge- (—e) spin- up holes created by ipi move in the same direction as the charge- (e) spin-up 
electrons created by ip].. The "constant" contains the truly constant ground-state energy of 
the spin down electrons, but also the term — f^e / 6{0)A{x) dx that subtracts a background 
electric charge. This charge gets discarded as we switch to the charge-conjugate picture in 
which charge-(— e) holes occupy the states that are not occupied by electrons. Keeping track 
of the "constant" restores the physical charge when needed. 

The vector potential A acts as a chemical potential and controls the location of the Fermi 
energy. In much of our discussion we will assume that when A = the Fermi energy lies 
midway between two edge-mode energy levels. This assumption is for illustrative purposes 
only. Indeed the detailed current /phase relation will depend sensitively on the exact location 
of the Fermi energy relative to the edge modes because varying 6 can make a level cross 
the Fermi energy, change its occupation, and cause a jump in the Josephson current. The 
sensitivity will manifest itself as Bohm-Aharonov oscillations in the Josephson current as a 
function of the magnetic flux through the Hall bar 

For our mid-spaced Ep we can make a gauge transformation to set A —t- at the expense 
of changing periodic boundary conditions to antiperiodic ones, and simultaneously redefining 
6{x). We assume that we have done this. The Bogoliubov-deGennes (BdG) equation for the 
eigenmodes is therefore 

-^Vd^ + \A{x)\e"'''^^^a^ ^ 




(4) 
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Equation (jl]) has a path-ordered exponential solution 

;;j^J)^e--Pexp{-./%.(o.«}(;;;°J). (5) 

where K{x) = |A is a hermitian matrix. Note that, in distinction to the 

usual BdG case, we did not double the number of degrees of freedom when we constructed 
the BdG operator, so all the BdG eigenmodes are needed. 

Only a part Q (the union of the two regions under the SC electrodes) of the perimeter of 
the Hall bar is in contact with the superconductor, and we set 

f/ = Pexp j-^y G SU(2). (6) 




As the perimeter of the Hall bar forms a closed loop, it was reasonable to impose periodic 
boundary conditions, but recall that these were changed to antiperiodic boundary conditions 
by the gauge transformation that removed A{x). The eigenmodes of the BdG operator 
Hamiltonian are therefore determined from the eigenvalues of U by requiring that 

= -e^^"^/-'*f/ r " I . (7) 

\VnJ 

Here L is the length of the Hall-bar perimeter. Now the eigenvalues of U will be of the 
form e^^'^ and so the energy eigenvalues are given by the requirement that (EnL/vd) ± = 
7r(2n + 1), or 

E„ = ^(7r(2n + 1)T0). (8) 

Note that if {u,v)'^ is an eigenvector of U with eigenvalue e**^ then —ia2{u*,v*) = {—v*,u*) 
is an eigenvector of U with eigenvalue e~^'^. Consequently if {un{x) , Vn{x))'^ is an eigen- 
function of the BdG operator corresponding to eigenvalue En, then (— f u* (x))"^ is an 
eigenfunction corresponding to energy —En- These facts follow from 

{ia2M-ta2) = -(T* =^ {ia2)U-\-ia2) = U, (9) 

and give rise to the usual antilinear S-wave BdG particle-hole symmetry "C" with = —Id. 
This symmetry must be distinguished from the approximate particle-hole symmetry arising 
from our linearization of the quantum Hall edge-mode spectrum. 

If the phase of the order parameter is constant in segments VLi^2 (the superconducting 
leads) then U = U2U1 where 

cos Da — ie*^" sin Da 

Here Da = \A\wa/vd where Wa is the width of lead a. The eigenvalues of U are e^**^, and, 
by taking the trace of U, we see that (p is given by the spherical cosine rule: 

cos (f) = cos Di cos D2 — cos 6 sin Di sin D2. (11) 

The spherical triangle (see figure [3]) arises because the matrices Ui and U2 are the spinor 
representations of successive SO (3) rotations through angles 2Di and 2D2 about axes sep- 
arated by the angle 6. It is shown in [l3| that such rotations can be combined through the 
use of mirrors that form the geodesic sides of the triangle. 
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FIG. 3: The spherical triangle that relates the eigen-phase (p to the order-parameter phase difference 
6 = 62 — Oi. 



From now on we understand by "0", the solution of (fTTj) that hes in the range < < vr, 
and by the vector (u, vY' the corresponding eigenvector of U . We similarly take "-En" to 
mean the combination 

En = j^ (2^(^ + V2) - 0) . (12) 
Now we make the Bogoliubov transformation 

In order not to over-count, we ensure that the modes are those that, after passing the 
superconductor, take the form (M„(a;), = e*^^"^"^'^^(M, f ), and (— f *(a;), m* (x)) = 

^-i{Enx/vi+<j>)(^_^* rjj^g Fermionic anticommutation relations coupled with the BdG 
eigenfunction completeness relations then require that 

{Ki, Ki) = {6„t, &mt} = {bni, bm^} = {b^^, b]^^} = 0, (14) 

and 

{bli, b^^} = {bl^, = Snm- (15) 
The Bogoliubov transformation simplifies H to 

oo 

H = ^ En{b\^bn^ - bnib\^) + const, (16) 

n=— oo 

the "constant" being the same one that was introduced earlier. It is not really a constant as 
it depends on the gauge field A, but it is independent of 6{x). Recall that the A dependence 
accounts for the total charge of the spin-down Fermi sea that was discarded when we made 



the particle- hole interchange for this spin component. The minimum-energy state is defined 
by the properties 



6„t|0) = 0, En> 0, 

6l^|0) = 0, En< 0, 

6„;|0) = 0, E^> 0, 

6ljO) = 0, E^< 0. 



(17) 



Using these, we compute the ground state energy to be 




(18) 



The quantity -Eground is formally divergent, but the physics resides entirely in the variation 
of -Eground with the phase difference 6 = 62 — Oi. Now as we vary 6 all En move in the same 
direction. The energy dependence on 6 largely cancels between the two sums. In order to 
extract the small, but non-zero, residuum we will have to regulate the sums in a controlled 
manner. This we do in the next section. 

III. COMPUTING THE CURRENT 

Given a Dirac-like spectrum of energy levels —00 < En < 00, the associated ground-state 
charge and current can often be expressed in terms of the spectral asymmetry [l^. This 
quantity is defined Q to be the regulated sum 



For energies of our form. En = Q;(27r(n + 1/2) — 0), a direct calculation shows that for 
— TT < < TT, we have 





(19) 
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FIG. 4: A plot of r]{(f)) showing the 27r periodicity. 

We may similarly define and compute an analytically-regulated version of the ground- 
state energy flTS]) : 



(£'ground)reg = Hm i - V Sgn(E„)E„e "l^"' + I 
I. n=— oo ^ 



02 vr 



a,---,, -7r<0<7r. 



This quantity also extends periodically outside the range — vr < (p < n — see figure El The 
subtraction needed for the existence of the limit is independent of and the constant — avr/G 
is the same as would be obtained by zeta-function regularization [17|. Let us also compute 



(ii^/gi-ound \ dcf 



JTo I - E «g-(^") (^) 

n=— oo ^ ^ 

foo 
a sgn(E„)e"''-^''l 



a- 



and observe that the regulated energy possesses the comforting property that 

^li^groundjreg - ^ ) ' 



(22) 



reg 



We will relate these energy derivatives to the ground-state expectation value of the diver- 
gence of the current operator. 
The current operator is 



3[x)= - 



SH 

6A{x) 



(23) 
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FIG. 5: A plot of a "^(-E'ground('?5*))reg showing the 27r periodicity 



If we include the contribution from the A dependent "constant" when taking the functional 
derivative, then the ground state current is 



= 2ev, ( K{x)\' + ^ \vn{x)A . (24) 

\E„<0 E„>0 ) 



If we ignore the "constant," the current becomes 

= eT;d(0|^|(x)V^t(x) - V';(x)V^j(x)|0) 

= efd ^ {\un{x)\^ - \vn{x)\^') - cvd ^ (|m„(x)|^ - |f„(x)|^) . (25) 

En<0 En>0 

These two currents differ only by the subtraction of + l"^"!^)!^) ^^e second 

case. This divergent sum is "5(0)" and independent of x by eigenvector completeness. When 
it comes to computing the current flowing in and out at the leads we can use either expression 
therefore. The second expression is the most convenient, and so we define 

(j(a;)),,g = lim | -ev, sgn(E„) (|u„(a;)p - |t;„(x)p) e"^!^"! I . (26) 

I. n=—oo J 

In our simple model |u„p(x) and |f„p(x) are independent of n, but do depend on whether 
X lies between the superconducting leads or not. This means that the edge-current differs 
in the two regions, and the difference is due to the Josephson current flowing in and out 
via the SC leads. We could compute \un\^ and |f„p in the two regions by diagonalizing the 
matrix U, but it is simpler, and more revealing, to relate the difference in the currents to 
the variation of the ground state energy with 6. 
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To do this we observe that 







\A\e~'^ -ivdid^ + ieA) 6+*^/^ 

-iVd{d^ - i{eA + x'/2)) |A|e^(^+'^) 

|A|e-*(^+^) ~ivd{d^ + i{eA + x'/2)) 

As the similarity transformation does not change the eigenvalues of the BdG operator, we 
see that 

E^[9,A] = E^[9 + X,eA + x'/2]. (27) 

The effect on the energy eigenvalue of changing 6{x) — t- 6{x) + 5d{x) is therefore identical 
to changing eA — )■ eA — (66)' /2. By first-order perturbation theory we compute the latter 
effect to give 



{n\SH\n) 

-Vd / dx [\Un{x)\^ - \Vn{x)\^) 5 A 



= J dxilUnix)]"^ - \Vn{x)\'^)-^S6{x) 

= -\vd I dx!^-^{\u4x)\'-\v4x)\')'^6e{x). 
Now, on combining this last result with equations (122|) and (125|) . we find that 



(28) 



ground / reg 



2e 



(V-j),,/^^(x)rfa; 



— -^Josephson(<^^2 ^ ^^l)- 



(29) 



Thus we see that the general result 

-^Joscphson ~ 

is consistent with our regularization scheme. 
From 



2c \ C^-£'ground 



h 



d9 



Vd ( 0' 



round 



L V27r 6 



we have 



d 



< (/) < vr. 



d 



-^Josephson 2c (-Zl/g]-omj(j)i.(;g 2e (^/gj-Qm^^j 



'reg 



dO 



(30) 



(31) 



(32) 



Figures |6|7|8] show how theses ingredients assemble to give the current/phase relation. 
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FIG. 6: A plot of the eigen-phase (p against 6 for the case Di = D2 = vr/2 — .2. We are enforcing 
the condition < < vr that is required by our Bogohubov transformation. 

d(/)/d0 



0.5 



-0.5 
-1.0 



1 




\ 1 2 3 

I - ^ — 


4 5 6 



FIG. 7: A plot ofd<p/de for Di = D2 = n/2 - .2 



To gain further insight, consider the case of "perfect couphng," where sin Da = 1 and 
= ±{02 — 6*1 + tt). In this case 



U 



-le 
162 g 



102 



—%e 

_Qi{02-ei) 















.g-i(e2-ei) 



(33) 



and so = (^2 — ^1) + tt. In the absence of relaxation, each 27r turn of 6 would put 
another particle into both the spin-up and spin down sea. In equilibrium however, the state 
ceases to be occupied as soon at its energy becomes positive. This change in occupation 
leads to a jump in the Josephson current as the state crosses the Fermi energy and its 
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■^Josephson (^)/(2e V^/L) 

0.6 




FIG. 8: A plot of /josephson/(2evd/L) against 9 for Di = D2 = vr/2 — .2. Observe how the 
discontinuities combine to give a smooth result. As D12 approach "perfect coupling" at Di = 
D2 = vr/2, the drops at 6 = 0, 27r steepen, and become level-crossing discontinuities. 

contribution is lost. The maximum possible current occurs just before or after the jump, 
and has /max = ±2efd/L. For ~ 10^ m/s and a perimeter of about 1 mm we get an upper 



bound on the 



and Zyuzin 



osephson current of about 1 nA. This is consistent with the estimate of Ma 



A physical picture for this upper bound is as follows: At the phase difference correspond- 
ing to the "jump," we have a spin- up/spin-down pair of levels lying exactly at the Fermi 
energy. At perfect coupling, the extreme equilibrium currents correspond to two possible 
cases: i) between the leads both zero-energy levels are empty whilst outside they are occu- 
pied, ii) between the leads both zero-energy levels are occupied and outside they are empty. 
Levels in the Dirac sea that are not at the Fermi energy cannot be left empty by a pas- 
sage under a lead, as this would lead to the energy being different in different regions and 
this is not possible in an energy eigenstate. Only the topmost energy level can contribute 
to the equilibrium Josephson current therefore, and this is the reason why the Josephson 
current is so small. To estimate its magnitude we note that in case (i), in each passage 
round the perimeter of the Hall bar, a pair of electrons passes from the Hall bar to the first 
lead and is returned to the Hall bar from the second lead. In case (ii) in each orbit a pair 
of electron passes from the first lead to the Hall bar, and is collected from the Hall bar at 
the second lead. This physical picture shows that the two possible Josephson currents are 
equal and opposite and have magnitude |/max| = 2efd/L. (Because it is easy to get confused 
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by Bogoliubov transformations, we provide, in Appendix A, a more detailed description of 
what happens to the particle content of the many-body eigenstates as they pass under the 
superconducting leads.) 



IV. COMPARISON WITH PERTURBATION THEORY 

The analytic regularization method used in the computations in the previous sections 
is standard in relativistic field theory but is perhaps less familiar in superconducting 
applications. As a check on its validity it is worthwhile (and non-trivial) to compare our 
all-orders in Di and D2 calculation with conventional perturbation theory. 

In the weak-coupling regime, where Di and D2 are small, the spherical cosine rule reduces 

to 

(j)'^ = Dl + Dl + 2D1D2 cose + 0{D^). (34) 
In this limit the ground-state energy and zero-temperature and Josephson current become 

^ground(^) = ^^{Df + Dj + 2D,D2 cosO), (35) 
L Zn 

and hence 

/josephson = -^^-C'iZ:'2 sin 6*. (36) 

We begin by verifying that fl35|) is correctly reproduced by the perturbation expansion. 

The Euclidean chiral propagator for zero temperature and anti-periodic spatial boundary 
conditions is 

wm.M,um = w^, - = J^-j-^J^ (37) 

where a, b =t, | and z = x + ivdT. The change in the A = ground-state energy due to the 
interaction 



H,^, = j \A{x)\ i^e''^^^^\{x)^lix) + e-^''^^V;(a;)^t(^)j dx (38) 
occurs at second order in lAL and is 



oo 



5i?ground = - dxj dx2 / rfr | A (xi) 1 1 A (xs) | e^^^'^^e'^^^-'^) (0 |TV'}(zi )^| (zi) ^^(^2)^t (^2) 1 0) . 

J J J —00 

(39) 

Here t = T2 — ti is the Euclidean time interval between Z2 and Zi. Now 

(0|T^{(zi)^|(zi)^;(z2)^t(^2)|0) = [G(zi - Z2)]' (40) 
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by Wick's theorem, and 
1 



4L2 



■ dr 



1 



(41) 



(sm[7r(a;i — X2 + ivdr)/L])^ "' \2T:Lvd, 
is independent of the separation x\ — X2 unless xi — X2 — ^ (mod L). The perturbation 

integral has four contributing regions: i) both xi and X2 in lead 1, ii) both x\ and X2 in lead 

2, iii) xi in lead 1, 2:2 in lead 2, iv) xi in lead 2, ^2 in lead 1. Recalling that Da — \/S.\wa/vd-i 

these combine to give 

1 



(^^^ground = vi{Di + + 2D1D2 COS ^) 



2TrLvd 



Vd 

2nL 



{Dl + DI + 2D1D2 cos I 



(42) 



This expression coincides with the weak coupling limit of the all-orders calculation. 

We can extend the comparison to non-zero temperature. At temperature T — f3~^, the 
Josephson current can be written as 

^Josephson - [j J (43) 

where F is the free energy. For a general spectral shift 0, we use standard methods to write 
down the partition function 



Z = exp{-/3F[0,/3]} 

/3vd 



= exp 



2n-l\2 



OO f 

J2 ^^p{" 



N=l 



27rL2 



wq 



(27rn + (py 



)^(l + «;-V"~')' 



(44) 



where q — exp{— TT^t'ci/L}, w — exp{— Pvd(t) / L} , and 

CXD 



n=l 



is the Dedekind eta function. We used the Jacobi triple-product formula to pass from the 
second line to the third. The sum in the expression for Z is squared because there are two 
independent Fermi seas (spin up and spin down) and their contributions to the partition 
function are symmetric under the interchange of (p with —0. By using the Poisson summation 
formula, we can rewrite the partition function as 



exp{-/?F[.#.,/)]} = (nm-^—^ 



E 



exp 



1 2ttL 



2vdf3 



n + inc 



VdP 
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[9s{<f>/2n\iL/vdP)y 



(45) 



^Josephson (T)//josephson (0) 
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FIG. 9: A plot of the effect of temperature on the perturbative Josephson current.. The horizontal 
axis is temperature in units of hv^/L. We see an effect as soon as the temperature becomes 
comparable with the 27r /itJ^/L level spacing of the edge energy states. 



Thus the free energy is given by 

F[(f),P] = c-^lnes{<P/27T\tL/v,f3). 
where c does not depend on 0. For small spectral shifts 0, we can Taylor expand 
F[0, /3] = c - ^0^^ In e,{(j)/27T\L/v,P) + 0(0^). 



(46) 



(47) 



We would now like to compare the expression f H7|) with that obtained by perturbation 
theory. At finite temperature the chiral propagator becomes 

1 e'{0\tVdP/L) es{z/L\tVdP/L) 



Here we are using the theta function definitions from jisj] , in which 



(48) 



ZT 



9,(z\r) 



J2 exp{ 

l= — oo 
oo 



+ 1/2)2 ^ 2vri(m + l/2){z + 1/2)} 
+ 2'Ki'mz^ . 



(49) 



m=— oo 



Thus 0{z\t) is odd under z -H- —z, while 9?j{z\t) is even. These properties were the ingre- 
dients used to assemble ( 148|) . which is specified uniquely by requiring the propagator to be 
analytic, doubly anti-periodic 



G{z + L) = -G{z), G{z + iVdl3) = -G{z), 



(50) 
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and for small z to obey 

G(.) ~ ^i. (51) 

It is this last property that makes it a Green function. 
In terms of Gi^z) we now have 

^^ground = - [dxi [ dx2 [^dT\A{x,)\\A{x2)\e'^'''h-''^''\Gixi-X2 + tVdT)f. (52) 



J J Jo 

The Xa integrals are the same as before, and, although it is little more complicated, the 
integral over r can still be evaluated in closed form. We begin by observing that [2TTiG{z)]'^ 
is analytic, has a double pole l/z"^ at the origin, is doubly periodic with periods ui = L and 
1^2 = ivdf^, and (from the 93{z\t) in the numerator) has a double zero at z = ^{coi + ^2)- 
These properties are sufficient to show that 

[2niG{z)f = p{z I wi, W2) - 63, (53) 
where p{z \ ijJi,ijJ2) is the Weierstrass elliptic function, and 

63 = p{{u}i + uj2}/2\uji,uji). (54) 
The Weierstrass zeta function is defined so that 

-^C{z\Ui,U2) = -piz\ui,uj2), (55) 
together with initial condition 

lim(c(^)--|=0. (56) 
We may therefore evaluate the r integral in terms of tabulated functions: 

/ [2niG{z)]'^dz = -((a + Ws) + C(a) - ^263 

J a 



-2ri2 -^263, 

■In 6^3(^1 - UJ1/UJ2) 



(57) 



z=0 



U2 dz"^ 

Here 27^2 = C('3' + W2) — Cl^^) = 2C(w2/2) is independent of a. The quantities in the second line 
of fl57j) are available in Mathematica"'"'^, and we use them to plot /josephson(7')//josephson(0) 
in Figure O 

It takes a little more work to obtain the logarithmic derivative appearing in the last line 
of ( 1571) . and so we relegate its derivation to Appendix B. Accepting that the claim is correct, 
and putting in the dimensionful constants, we confirm that our all-orders evaluation of the 
free energy coincides with the perturbation theory calculation in the weak coupling regime . 
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V. DISCUSSION 



We have shown that the maximum possible Josephson current for a pair of spin- up/spin- 
down QH edge states is rather small for typical Hall bar geometries. The bound is small 
because the relevant length and energy scales are set by the perimeter of the Hall device 
rather than the separation of the superconducting probes. Also, unlike a typical Josephson 
device, there is only one conduction channel per pair of edge modes. This last observation 
means that nothing is to be gained by making the superconducting leads overlay deeper into 
the Hall bar. 

It may seem strange that we have so far discussed quantum Hall physics with no mention 
of the magnetic field that is necessary for its existence. The field, however, has only a few 
consequences for our discussion. Obviously the superconducting leads must be constructed 
of materials that remain superconducting in a field of few Tesla at temperatures of about 
IK, but this is not hard to achieve. The leads must also be narrow enough that the order- 
parameter phase does not vary widely within the part of the lead that is actively coupled 
to the 2DEG. A subtle point in this regard affects the claim that the Josephson current is 
independent of the separation of the leads. The phase difference that we have equated 
to 62 — Oi should be understood as the gauge invariant quantity 6 = 62 — Oi — 2e J^^ A dx. 
Now a quantum of magnetic flux lies between each of the edge-state energy levels and if the 
effective is not to vary with the energy level index n, only a small fraction of this flux 
should pass between the leads. The leads should not be spaced apart by more than a small 
fraction of the perimeter. 

An effect that we have not considered here, and one that may well allow for larger 



currents, is "edge reconstruction" 



20 



221 ]. A reconstructed edge, with its alternating strips 



of compressible and incompressible 2DEG can allow many more levels to lie exactly at 
the fermi energy and so have their occupation number changed without a change in energy. 
These levels have zero drift velocity, however, so it unlikely that they contribute significantly 
to the Josephson current. 
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VII. APPENDIX A 



The maximum possible Josephson current occurs when we have both perfect coupling 
(sinDi = sin 1^2 = 1) and cos 6* = 1. In this special case we have 



U,=U2 



-i 
-i 



U = U2U1 



2ninx/L _|_ 



The Bogoliubov mode-expansion f[T3|) then becomes 
for X in region (1), and 

"I 



1 
1 



1 /o 



(58) 



-l-KinxjL 



\fL \ 1 



(59) 



i\)\{x) 



2-iTinx/L _|_ r t 



Vl \ 



(60) 



for X in region (2). (The numbering of the regions refers to figure [H) 

In these mode-expansions, the operators 6„-f- and annihilate or create quasiparticles 
with energy \En\ = 27[Vd\n\/L. We compare these expansions with the free-particle plane 
wave expansion 



tfj^ix 



^2Trinx/L _|_ 







'2ninx/L 



Vl \ 1 



(61) 



where the operators a„-|- and a}^^ annihilate and create electrons. We see that we can identify 



Ki = cini, = all 



(62) 
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in region (1), and 

but = ialnp ^It ^ 

Ki = -iai^, = +ia^nt (63) 

in region (2). We now use these identifications to examine what happens to the particle 
content of the many-body eigenstates as they drift under the superconducting leads. 

We first note a minimum-energy eigenstate must be annihilated by bn-^ and bni for n > 0, 

and by bl^^ and 6^^^ for n < 0. Let us define the eigenstate |0) by requiring that it is killed 
by all these operators, and also by 604, and 6of. Then the states 

|0), bl^\0), 6;jO), bl^bl^lO), (64) 

all have the same energy, making the ground state four-fold degenerate. 
With the operator identifications established above, we find that 

-1 

|o)= n («iA)N^pty) (65) 

n=— 00 

when X lies in region (1), but in region (2), where 6ot ^04, are identified with a^^ and a^^ 
respectively, we must have 

-1 

io)««Wt n («Wt)i®"^p^y)' (66) 

n=— 00 

for it still to be annihilated by 6ot and 604,- We see that the occupation number of the 
energy levels for n < are unchanged, but |0) picks up a pair of n = electrons from the 
superconducting lead as it passes under it. Similarly the state feJ^^feJ^jO) loses a pair from the 
n = level. 

The state bl^\Q) is annihilated by and oq^ in region (1), and these become respectively 
004. and in region (2). The particle content of this state is unaffected by its passage under 
the lead therefore. Similarly bl^\0) retains its particle content. 

Now consider an excited state, for example 6]^-|-6|j-|-|0) with m > 0. This state has energy 
E — 2'KVdm/L. In region (1) it has particle content 

-1 

n («l4«lt)N"^pty)' (67) 

n=— 00 
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and so consists of a Dirac sea together with an electron in a positive energy level. In region 
(2) it becomes 

-1 

n=— oo 

which consists of a Dirac sea which has lost an electron from a negative energy level. After 
passing the superconductor therefore, the state has the same energy and spin, but the 
electron has become a hole. 



VIII. APPENDIX B 



We wish to establish the third line of (157|) . which reads 



/ 

•J a 



a+u)2 ]^ ^2 

{p{z\uJi,UJ2) - es)dz = — — ln6'3(2;| - wi/wa) 



U2 dz^ 

This result follows indirectly from the related integral 



(69) 



2 = 



J a 



{p{z\ui,uj2) - e'i}dz = -2?7i-wie3 

1 9'^(0\uj2/(-Jl] 



ioi 6'3(0|a;2M) ' 

\n93{z\uj2/uJi) 



. (70) 

Ui dz^ 

Here we require Im(u;2/wi) > for the theta functions to converge. To establish (ffOl) we 
observe that second line follows from the first by combining two standard formulae: 



18| Eq 5.2), and 



I 3 ^'fO|r) 03(O|r) 
1 l^"'(0|r) 



2r/i = - 



Wi3 ^'(0|r) ' 



(71) 



(72) 



19| §21.43.) Here r = UJ2/UJ1 with Imr > 0. The third line of (!70|) follows from the second 



because e^(0|r) = 0. 

To derive ( |69l) however, we need the integral over the uj2 = ivdP imaginary period, and 
not over the ui = L real period. Because of the positivity condition on the imaginary part 
of r, we cannot change the integration path by merely interchanging Ui ^ U2 in equation 
rOj) . We need to be more subtle. By changing (c<;i,c<;2) — )■ (— (X'2,(X'i) in ( [70|) . we obtain 



l_d^ 
UJ2 dz'^ 



ln6'3(2;| - U1/U2) 



2=0 



a—W2 



{p{z\ -U2,uji) -esjdz. 



(73) 
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This last equation is legitimate because Im(c<;2/i^i) > implies that lm{—ijji/u2) > 0. We 
now manipulate 

RHS = - / {p{z\ - uj2.,uJi) - e-i} dz 

J a 

= -/ {p{z\uji,uj2) - 63] dz, (74) 

J a 

where the last line follows from the invariance of p{z\ Ui, U2) under modular transformations 



uj[\ fa b\ / UJi\ fa h 

uj'2 ) \c d ) \UJ2 ) \c d 



gSL(2,Z). (75) 



From this we immediately deduce eq. 
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